function jv = Jvec_TDIP(invP,etaModel,j,v)


q = Geta_vec(invP,j,v);
w = calcAhInv(invP,etaModel,q);
jv = -projHs(invP,w);

